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A solenoid oscillating in vacuum will pair produce charged particles due to the Aharonov-Bohm (AB) 
interaction. We calculate the radiation pattern and power emitted for charged scalar particles. We 
extend the solenoid analysis to cosmic strings, and find enhanced radiation from cusps and kinks on 
loops. We argue by analogy with the electromagnetic AB interaction that cosmic strings should emit 
photons due to the gravitational AB interaction of fields in the conical spacetime of a cosmic string. 
We calculate the emission from a kink and find that it is of similar order as emission from a cusp, 
but kinks are vastly more numerous than cusps and may provide a more interesting observational 
signature. 



Fifty years ago Aharonov-Bohm (AB) [T[ showed that 
charged particles can scatter non-trivially in the pure 
gauge potential outside a thin solenoid. The effect has 
been experimentally investigated in various setu ps B, 3 
and quite recently (theoretically) in graphene [f 3I |. In 
cosmology, AB scattering is relevant to the dynamics of 
cosmic strings, as they move in the ambient cosmologi- 
cal medium [5|, [f| and the AB interaction could be a link 
between dark strings and the visible sector Q. 

AB scattering is purely quantum, having no classical 
analog, and also purely topological, as it arises from the 
non-trivial change in the phase of the wavcfunction of 
a charged particle as it is taken in a closed path around 
the solenoid. In this paper we study if the AB interaction 
can also lead to particle creation - if a thin long solenoid 
oscillates, does it produce electron-positron pairs? Such 
radiation would be purely quantum and purely topolog- 
ical. 

The question of AB radiation has been partially ad- 
dressed by Alford and Wilczek Q where the problem is 
set up (but not solved). In this paper, we shall take up 
this task using two approaches. The first is to use the 
AB phase, defined as e = e<!> where e is the charge and <& 
the magnetic flux through the solenoid, as a small con- 
trol parameter. Here the calculation can be done using 
the Feynman diagram language with the gauge potential 
of the moving solenoid treated as a classical background. 
The method can also be applied to find the AB radiation 
rate from oscillating cosmic string loops. The disadvan- 
tage of this perturbative approach is that we lose track 
of the expected periodicity of all quantities in the AB 
phase. For example, the classic "sin(7rev)" behavior of 
AB scattering, where a = e/27r, does not appear in this 
treatment. Our second approach is to allow the AB phase 
to take any value but to only consider slowly oscillating 
solenoids i. e. we expand in the velocity of the solenoid. 
The advantage now is that we explicitly see the periodic 
dependence of the radiation rate on the AB phase. How- 
ever, we cannot treat relativistic motion such as that of 
cosmic string loops in this approximation. We can, of 
course, take both the AB phase and the oscillation speed 
to be small, in which case both methods should agree, 
and we show that they do. 



Another novel aspect of AB radiation is that it may 
also apply in the gravitational context. In a conical met- 
ric, such as that of a cosmic string, the wavefunction of 
any particle that goes around the string acquires a phase 
equal to the conical deficit, which we denote S = 8nGfi 
where G is Newton's gravitational constant and /x is the 
string tension. Then, if a string oscillates, it emits all par- 
ticles including, for example, photons. This may make 
cosmic strings, even those whose constituent fields have 
no interactions with photons, visible in the electromag- 
netic domain. We will discuss gravitational AB radia- 
tion in Sec. |V] where we rely on the analysis developed 
in the earlier sections. Although we arrived at gravita- 
tional AB radiation via the gauge AB process, the effect 
is just that of particle production in the time-dependent 
gravitational background of a cosmic string considered 
by Garriga, Harari and Verdaguer in Q (also see @). 



The rest of the paper is organized so that we first dis- 
cuss the gauge potential produced by a moving solenoid 
in Sec. HI as given in Q. Then we find the AB radiation 
from an infinite, straight, oscillating solenoid in Sec. IIII1 
using first the small AB phase approximation and then 
in the slow velocity approximation. In Sec. IIVI we cal- 
culate AB radiation from cosmic string loops, providing 
two explicit examples, one of a loop with kinks and the 
other of a loop with cusps but no kinks. In Sec. [V] we dis- 
cuss gravitational AB radiation. We conclude in Sec.lVIl 
The formalism to find AB radiation in slowly changing 
backgrounds is set up in Appendix [3] 



In this paper we have restricted ourselves to AB radi- 
ation of spin zero (scalar) particles and are mostly con- 
cerned with the case when the emitted particles are mass- 
less. The extension to fermions is conceptually similar 
but the spin orientation provides another degree of free- 
dom and makes the calculations more involved at a tech- 
nical level. We will address the fermion calculation in a 
separate publication fic| . 
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I. GAUGE POTENTIAL OF MOVING 
SOLENOID 

The first step in calculating AB radiation is to find 
the gauge potential for a moving solenoid. This step was 
taken in Alford and Wilczek Q. 

A moving solenoid with magnetic flux <f> will be a cur- 
rent source in Maxwell's equation, and the current can 
be written as 



t($) _ £ 



(1) 



where 



S af) {x) = J dTda^e ab d a X a d b X fs 8 ii) {x- A>,t)) 
drda{X a X ' - X p X a ')S {4) (x - X(a, r)) (2) 



where a,r are world-sheet coordinates, X fJ '(a) denotes 
the position of the string, overdots and primes denote 
derivatives with respect to r and a respectively, and "f a b 
is the worldsheet metric, with 7 denoting its determinant. 
In what follows, we shall use r = t, while a = z for 
a solenoid along the z— axis. In applications to cosmic 
strings, we will also use the gauge conditions 

X ■ X' = , X 2 + X' 2 = (3) 

The field strength of a static, thin solenoid satisfies 

d»F^ = J v (4) 

and the solution is 



rp ^ , qafi 

~ 2 t M zyQ P' J 



(5) 



The corresponding gauge potential can be written in 
Lorcnz gauge 



$ 1 a 
A — _ e „f)^ Q a P 



(0) 



where d 2 is the D'Alcmbcrtian operator. We will only 
need the gauge potential in momentum space, 



"i" — '2 t M I/a P ^,2 



(7) 



where overtilde's denote the Fourier transformed variable 
S Q/3 (fc) = / d 4 x e^^^V) • (8) 



Note from Eq. ([5]) that the field strength vanishes ev- 
erywhere except at the location of the solenoid and so 
a moving solenoid does not radiate classical electromag- 
netic waves. 
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FIG. 1: Feynman diagram for pair creation by oscillating 
solenoid. The solenoid creates an oscillating gauge poten- 
tial, A^, while the field strength is always zero outside the 
solenoid. The time dependent pure gauge leads to parti- 
cle production. 



II. AB RADIATION 

Consider the interaction term J^A^ where, for exam- 
ple, for a complex scalar field, \i we have 

Jfi = ie(x*dftX ~ xd^x*) (9) 

where e is the charge of \- The lowest order amplitude 
for x production is found from the diagram in Fig. [1] and 
can be written as 



(f> 

iM = -i- efiua ^S^(k)J v {p,p')\ k=p+pl (10) 



with J u (p,p') = ie(p-p r y. Then 

... e (p + p'Tip-p'Y 
lM - 2 W (p+pO 2 



« — yP F a 'V"'(P+P')'^ 



x / drda{X a X p - X p X°")e 



(11) 



where e = e$ is the AB phase and is assumed to be small. 

The number of particles produced is given by (see, for 
example, Eq. (4.74) of [n|) 



dN = 



d 3 p 1 dV 1 , 

-\M\ 2 (12) 



(2tt) 3 2w (2tt) 3 2u' 



where uj = \/p 2 + m 2 (similarly ui') and m is the mass 
of the scalar particle. 



III. AB RADIATION FROM OSCILLATING 
SOLENOID 

In this section we will consider AB radiation from 
an infinite straight solenoid that is oscillating back and 
forth. We will do the calculation in the small AB phase 
approximation. The oscillating solenoid is also treated 
in the non-relativistic approximation in subsection IIIIB1 
and the results of the two methods agree in the regime 
where both are valid. 

We will take the position of the solenoid to be described 
by 



X 1 -' = (t,Asm(nt),0,a) 



(13) 
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where A is the amplitude and fl the frequency of oscilla- 
tion. Note that X M is not the solution to any equation 
of motion but just a convenient oscillatory choice that 
provides us with a background in which to study AB ra- 
diation. 



A. Small AB Phase Approximation 



J{n, x) 



Inserting Eq. (fT3|) in Eq. ([TTj) gives 



2e 



iM = ~ (p + p') 2 ^ PxP ' y ~ PvP '^ Iq + ^' Pv ~ up ' y)I 



where w = p°, uj' = p ' 

1° = I drda e -^+P'> A " 



ij 

(14) 



I 1 = J drda QA cos(f2r)e 



-i(p+p')-X 



(15) 



(16) 



The integration along the straight string is trivial, and 
leads to a 2ttS(p z + p' z ) factor in 1° and I 1 . 
To simplify further, we use the relation 



(17) 



which can be checked using the integral representation 
for the Bessel functions 

J„(a) = l- e -i(nr- aS m T ) dT (lg) 

Then 

/0 = E UP 1 A)5{n (P°/n))S( Pz + p' z ) (19) 



and 



I 1 = 2v 2 AY J [Jn-i(P l A) + J n+1 (P l A)] 

n 

x6(n-(P°/n))6(p z +p' z ) 



(2n) 



{- J2 nJ n (P 1 A)S(n - (P Q /il))S( Pz + p' z ) (20) 



P 



where P = p + p 1 and J n (x) is the Bessel function of n th 
order. 

We can now insert the expressions for the integrals 1° 
and I 1 in Eq. (JUJl to get 



iM = -8n z e 



2 (Py-Py) 



oo 

x Jn((p x +p' x )A)6(u; + u>' - nn)5{ Pz + p' z ){2l) 




FIG. 2: Plot of J n (x)/x for n = 1, 2, 3, 4, where the amplitude 
of the peak decreases with increasing n. 



At first glance the amplitude depends on six variables, 
namely the six different momentum components. How- 
ever, Eq. (f2"T|) shows that it only depends on the combi- 
nations p x +p' x , p y —p'y and p z +p' z as long as the overall 
energy is a harmonic of the oscillation frequency. Note 
that the amplitude is not singular at p x +p' x = because 
the Bessel functions are proportional to (p x + p' x ) n at 
small argument and the sum in Eq. (|2~T|) starts at n = 1 
(see Fig.[5]for plots of J n {x)/x). The amplitude is largest 
when p y — p' y is maximum and p x + p' x = = p z + p' z . 
Together with the total energy constraint 



uj + uj' = \Jp 2 + m 2 + y p' 2 + m 2 = nft, (22) 

we find that the maximum amplitude is along the y di- 
rection, with p x = = p' x . 

The number of particles emitted follows from Eq. (|12l) . 
On squaring Ai we get a double sum but the energy con- 
servation delta function reduces the expression to a single 
sum and introduces a 6(0) on the right-hand side. Divid- 
ing out by 27r<5(0) gives the rate of particle production, 
denoted by N, 



d 3 p 1 d 3 p' 1 (87r 2 e) 2 fp y - p[ 



dN 

(2tt) 3 2uj (2tt) 3 2uj' 2tt \p x + p' x , 

OO 

x I J "((^ + p' x )A)\ 2 5{u + u' - nn){5(p z +p' z )} 2 (23) 

The integration over p' z leads to another factor of (5(0). 
Dividing out by another factor of 27r<5(0) gives the rate of 
particle emission per unit length of the solenoid, denoted 
byN', 



dN' 



16tt 4 



dp z d 2 p±d 2 p'j 



, 1 fPv -Py 



ojuj' \p x +p' x 



Y,\ J n{iPx+p' x )A)\ 2 + (24) 



where the _L subscript denotes the x, y components. Also, 
now uj' = yp'j 2 + P 2 + rn 2 . 
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The expression for the rate of particle emission depends 
on five of the momenta. The dependence on p z enters via 
lo and uj 1 , and it is clear that the largest emission is for 
p z = i.e. in the plane perpendicular to the solenoid. 
So we restrict our attention to p z =0 and also to the 
massless case, m = 0. These assumptions imply that 
lu = \p± | , lj' = \p', I , and we obtain in polar coordinates 



dN' 



(Pz = 0) 



^—^d P± dp' ± d9d9' ' v Py 



dp z 16tt 4 \Px+Px, 

OO 

x \Jn(( Px +p' x )A)\ 2 S( P± +p' x - nil) (25) 

n=l 

where p x = p± cos9 etc.. 

The integration over p' ± can be done giving the emis- 
sion rate per unit length in the orthogonal plane as a 
function of the three variables, p±, 9 and 9'. The angu- 
lar distribution is of greatest interest and is shown for 
n = 1 and n = 2 in Fig. [3l taking p± = nil/2. The figure 
shows back to back emission, with the n = 1 emission 
being dipolar and maximum along the y direction. 

The total radiation rate can be found by integrating 
Eq. (|24p over momenta and also doing the sum over har- 
monics. The sum is convergent because the Bessel func- 
tions at fixed argument fall off exponentially with n for 
large n. We will work only in the to = case and restrict 
attention to particle emission in the orthogonal plane, as 
in Eq. |25|) . This is equivalent to working out the prob- 
lem in two spatial dimensions, and we shall denote the 
particle emission rate by N^ 2d K The delta function can 
be integrated out by going to polar coordinates. We get 



7\T(2d) 



e 2 Q 
16tt 4 



^ nC n (a n ) 



(26) 



where a n = nflA and 
-l 



C n (a) 



dq 



o 



d6 f' d0> (^-(l-gjBinf 
-it J -7i \qcos9 + (1 - q)cos9' 



'\\|2 



x\J n (a(q cos 9 + (1 — q) cos0'))| 



(27) 



These integrations have to be done numerically. 
The energy radiation rate in the n th harmonic is 

E n = nVLN n 



(28) 



and Fig. 0] shows a plot of ln(£„) versus n, where £ n = 
E n l6ir 4 /(e 2 n 2 ) for ClA = 0.1, 0.5, 1. Note that the energy 
falls off rapidly with increasing n and so essentially all of 
the particles are emitted in the n = 1 harmonic i. e. with 
energy f2. 

The total energy emission rate £ as a function of the 
parameter flA is shown in Fig. El 

We have not calculated the emission rate for massive 
particles, to 7^ 0, in detail. From Eq. (f2"5]l we see that 
to enters through uj and lo' in the phase space factor 
and then again in the energy conservation delta function 




FIG. 3: Here we show the radiated power as a function of 8 
and 9'. The top figure is for n = 1 with \p\ — nQ/2 — \p'\ 
and the lower figure is for n — 2. There is maximum radiation 
along the line 9 + 9' = n, implying back to back emission. For 
n — 1, the peaks occur at 9 — n/2, 3n/2. For n — 2, they 
occur at 9 = tt/4, 3ir/4, 5tt/4, 7tt/4. 



but nowhere else. The delta function can be non-zero 
only if nQ, > 2m i.e. n > 2m /CI. If the oscillation fre- 
quency is small compared to the mass of the particle - 
as is relevant for oscillating solenoids in the laboratory 
when the charged particle is the electron - the Bessel 
functions decay exponentially with large n and we ex- 
pect the emission to be exponentially suppressed. The 
exponential suppression is also seen in the plot of Fig. [4] 
It is noteworthy that the total energy of the pairs pro- 
duced add up to nfl by virtue of the delta function con- 
straint in eq (24). Commonly in quantum mechanics, 
perturbations at frequency produce excitations with 
energy hCl (as for example in optical transitions in atomic 
physics). Here excitations with energies equal to all mul- 
tiples of the fundamental frequency are produced (al- 
though in practice the higher harmonic production may 
be small except in an ultra- rclativistic limit). This sit- 
uation is also encountered in laboratory scale phenom- 
ena such as Mossbauer spectroscopy and the excitation 
of conduction electrons by the motion of impurities in 
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B. Small Oscillation Frequency Approximation 



FIG. 4: Ln(£„) vs. n for Q.A = 0.1 (lowest curve), 0.5 (middle 
curve) and 1.0 (top curve). 
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FIG. 5: £ vs. QA for the oscillating solenoid. 



metals. Mathematically the appearance of excitations at 
higher harmonics is traceable to the Fourier expansion eq 
(17). 

In the next sub-section, using the formalism developed 
in Appendix \^ we re-analyze AB radiation from an os- 
cillating solenoid for arbitrary values of the AB phase, e, 
but in the limit of small HA. We recover the periodicity 
of the result as a function of AB phase. We also show 
that the results of the two different methods agree in the 
limit that both e and QA arc small. 



Now let us consider the circumstance that the flux is 
not a small perturbation. In this case we can still obtain 
analytic results for particle production in the limit that 
the solenoid motion is non-rclativistic. These results are 
non-perturbative in the flux and show the flux periodicity 
expected of the Aharonov-Bohm effect. Moreover, the 
small flux limit of the results obtained in this section 
match the low velocity limit of the results in the previous 
section. 

The essential physics here is that for a stationary 
solenoid the eigenmodes of the scalar field depend (peri- 
odically) on the flux $ of the solenoid. If the solenoid is 
moved to a different location the eigenmodes and hence 
the quantum vacuum is altered. If the solenoid is moved 
slowly from the initial to the final configuration, in the 
adiabatic approximation it will simply pass from the ini- 
tial vacuum to the final vacuum, but post-adiabatic cor- 
rections lead to particle production. The technical imple- 
mentation of this idea is relegated to Appendix [X] where 
we analyze the leading post-adiabatic dynamics of an as- 
sembly of oscillators with slowly varying eigenmodes. 

We will work in cylindrical co-ordinates and in the 
"singular gauge" where the gauge potential of a static 
solenoid is 



A" = (o,o,$e(x)<%),o) 



(29) 



Here O is the Heaviside function. The gauge poten- 
tial vanishes everywhere on the xy— plane and is singular 
along the positive x— axis. It is easy to check that the 
line integral of the gauge potential along a path enclosing 
the origin yields the flux $. 

The eigenmodes, ipkqi, for a scalar field in the back- 
ground of a solenoid along the z— axis have been found in 
Ref. pj. We will use those eigenmodes in singular gauge. 
In order to box normalize the eigenmodes, we adopt pe- 
riodic boundary conditions along the z-axis and take the 
length of the solenoid to be L. Along the radial direc- 
tion we find it convenient to impose Dirichlct boundary 
conditions at r = a (i.e. ipkqi(a) = 0) if £ is even, and 
Neumann boundary conditions (i.e. d r ipkqi(a) = 0) if I is 
odd. These boundary conditions are convenient because 
they allow us to do certain integrals in closed form; they 
also give the physical result that the particle production 
vanishes in the limit that the solenoid flux goes to zero. 
The values of k now take on discrete values for finite a, 
but the set becomes a continuum in the limit a — >• oo. 
Then the singular gauge eigenmodes are 



(30) 



where < 9 < 2ir, and a = e/2ir. The modes are labelled 
by the radial momentum k, the z-momentum, q, and the 
azimuthal quantum number I. The frequency of the mode 
is 



(31) 



G 



Note that the eigcnmodes have a discontinuity along the 
positive x-axis. This is because we are working in sin- 
gular gauge where the gauge potential is singular on the 
positive x— axis. 

The normalization integral of the modes is 



r 2ir 

dr r I d9 



dz ^k'q't'^kql = 5ll>5qq>5kk<- (32) 



where we have taken ka^> \(l — a) 2 — 1/4| and used the 
asymptotic form of the Bessel functions 



Jv{x) f 

together with the integral 
dr r J v {kr)J v {k' r) 



vix 

T 



(33) 



— [(Mka)) 2 - J^_ 1 (A;a)J I/+1 (fca)]5 fc fc' (34) 

As in the previous section we assume that the solenoid 
oscillates in the xz— plane while remaining aligned with 
the z-axis. The modes for a shifted solenoid arc obtained 
by appropriately translating the modes eq (|30[) . Hence 
the transition element designated ap a in the appendix, 
Eqs. (|A6|) and (|A15|) . is here given by 



/•a /■27T 

ak'q'i\kqi = drr d9 
Jo Jo 



L d 

dz ^'i'Vq^W- 



A straightforward calculation reveals 

2k sin(7ra) / k' 



O-k'q'l'.kql — —L5 q 

for 1 = 0,1' = 1 and 

0>k'q'l',kql = — L5, 



k' 2 - k 2 



2k sin(7ra) / k' 
qq ' k' 2 -k 2 



(35) 



(36) 



(37) 



for I = 1, I' = 0; for all other combinations of I and I' it 
vanishes. 

Using Eq. (|A16[) and taking the continuum limit 
L, a — > oo we obtain dN' = dk dq I(k, q) where 



I{k,q) 



v 2 sin 2 (not) k k 2 



16tt 2 



LJkqQ 2 



2o 



2a-2 



x<d(tt - uj kq - \J q 2 + m 2 ) 



(38) 



is the rate of particle production at wave- vector (k,q). 
(To include antiparticlcs, we would double this rate.) 
Here k c = (fl 2 + k 2 - 2flu} kq ) 1 / 2 . Eq. applies for 
< a < 1; by virtue of the periodicity of the scalar field 
modes, I(k,q) is a periodic function of a with period 1. 

Some of the quantities in Eq. (|38|) have a straightfor- 
ward interpretation in terms of the kinematics of pair 
production. Let the pair of particles have momenta (fc, q) 



and {k 1 ,q'). By conservation of z-axis momentum q+q' = 
0. By energy conservation Wfc g = Q— LJk'q'- The minimum 

value of ujk'q' is \/ q 2 + m 2 ; hence uik q < — \J q 2 + m 2 ; 
this accounts for the Q function in Eq. Similarly if 

one particle has momentum (fc, q), then k c represents the 
radial momentum of the second particle. 
In the limit a — >• Eq. ([55)1 simplifies to 



16 \ V 0J k q 

(39) 

To make contact with the results of the previous sub- 
section we consider the non-relativistic limit of the re- 
sult for dN', (Eq. (|24p) in which only the first harmonic 
contribution is significant. If we write p in cylindrical 
co-ordinates as p = (k,9,q), integrate dN' over p' and 
9, halve the result since the relativistic expression counts 
production of particles and antiparticles, we obtain the 
flux in Eq. (f5!)|) . Thus the non-relativistic limit of the 
small flux approximation coincides with the small flux 
limit of the non-relativistic moving solenoid approxima- 
tion, providing a useful check on both calculations. 



IV. AB RADIATION FROM COSMIC STRING 
LOOPS 

A cosmic string loop trajectory is written in terms of 
left- and right- movers 



X = \[a{a-T) + b{a + r)] 

The choice of gauge, Eq. ([3]), requires 
\a'\ = 1 = l&'l 



(40) 



(41) 



where primes denote derivatives with respect to the ar- 
gument. 

Insertion of Eq. (j40|) in (fTTj) with some simplifications 
gives 

e (p+p'Yip-p'y 



where 



i a 



da+b a 'e- lk - b/2 
da^a a 'er lk - a ' 2 



(43) 



and a± = a±r and k = p+p' . Note that the integrations 
have been extended to the entire a, t plane. This is valid 
since the integrands are periodic functions in both a and 
t, and the only effect of the larger domain of integration 
will be to yield Dirac delta functions instead of Kronccker 
delta's. 

The integrals, I± satisfy the orthogonality relations 



k ■ I± = 



(44) 
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which allow us to relate I± to I±, 

k I ± . 

By using these relations, we find 



k°I° ± 



J — p p J + — 

kitk 



2fc° 



■(p-p')-(l + xl_) 



Now we can write 



(45) 



(46) 




FIG. 6: Three snapshots of a degenerate kinky loop. 



A. Kinky Loop 



4(^)2 



1ST 



16(w + to') 



— \(p-p').l + xl_ 



(47) 



The integrations in Eqs. (|43| extend over an infinite 
range whereas the loop dynamics is periodic. We can use 
the series representation of the Dirac delta function 



-, +oo 



-7T < X < +7T 



(48) 



to write the integrals over only one periodic domain 

oo 

il = £ 



n=l 



V 4tt 



= J2 ^ ' dcT-a a 'e- ik - a ' 2 s(^- 



n=l 
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The integrals 7± arc non-vanishing only for discrete pos- 
itive values of the total energy 



U) + U) 



Ann 



(50) 



where n is a positive integer. Therefore it is convenient 
to write the squared amplitude for a fixed harmonic 



\M\; 



167m / 



V 



x/ («)|2 



(51) 



where 1^ refer to the integrals in Eq. (|49)) together with 
the constraint in Eq. (|50"j) . 

The squared amplitude can now be written as 

|A4| 2 = yf;|M|^(^-nWo) (52) 



where <5_e(0) is an energy conservation delta function and 
we will divide out by 2it6e(0) to obtain particle radiation 
per unit time. 

We will now find the radiation in the m = case 
from two different kinds of loops, one that contains kinks 
(sharp corners where the tangent vector is discontinuous) 
and the other that contains cusps (points on the loop that 
momentarily reach the speed of light). 



Here we consider the "degenerate kinky loop" 

a = a_A, < cr_ < L/2 

a = (L-a-)A, L/2<a_<L 

b = <t + B, < er+ < L/2 

b = (L-a+)B, L/2<cr+<L (53) 

where A and B are two fixed unit vectors. Also note 
that 



6° 



(54) 



so that X° = t. 

This loop is "degenerate" because it consists of four 
straight segments and is "kinky" because of its four cor- 
ners. The four straight segments propagate with constant 
speed but shrink and expand due to the motion of the 
kinks (see Fig. [B]). 

Now the relevant integrals can be evaluated in the m = 
case and with k° = Airn/L, 



r(») 



- AY Bin((l + *-A) T )A (55) 

J + = =Tj-sin( l-K'B — IB 56 

nir 1 — [k ■ B) z V 2 / 

where K = k/k°~ L(p + p')/Aim. This gives 

16 V 



(p-p')-Ax B\ 2 

sin 2 ((1 + k ■ A)*f) sin 2 ((1 - k ■ B)*f) 
[1 - {k- Af] 2 [l - {k-B) 2 ] 2 " 



(57) 



The expression in (|57l) is finite as can be seen by taking 
limits, say as k ■ A — > 1. This shows that kinks on loops 
give off finite AB radiation at every harmonic. (The sum 
over harmonics, however, will be seen to diverge.) Maxi- 
mum emission occurs when p—p' is parallel to Ax B i.e. 
perpendicular to the plane of the loop. Also, the factors 
depending on k are maximum for k- A = = k-B, imply- 
ing that k, hence p + p' , is perpendicular to the plane of 
the loop. Since the amplitude vanishes for p p' = 0, we 
conclude that the maximum emission occurs for p = —p' 
and perpendicular to the plane of the loop. 
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The differential rate of particle production is given by 
Eq. CCU) 



dN = 



1 d 3 p 1 d 3 p' 1 
2L (2tt) 3 2^ (2tt) 3 2^7 



x V|AC=o ; n<K«-(^ + ^)L/47r) (58) 



n=l 



The total energy emitted is 



oc 



(59) 



where 



• 2n f d 3 p 1 f d 3 p' 1 

" ~ L2 n 7 (>p2W (2^277 

x <5(n-( W + a;')i/4vr) 



(60) 



The integral can be evaluated by numerical means. 

However it is prudent to check the scaling with n be- 
fore performing detailed calculation. Due to the energy 
conservation delta function we estimate u> ~ n, uj' ~ n, 
d 3 p ~ n 3 , c? 3 p' ~ n 3 , and the delta function contributes 



The contribution of the factor \AA\ 



i 

m— 0;n 



is less 



clear: a factor of n~ 4 is left over after re-scaling all the 
momenta and there may also be factors of n due to con- 
tributions from the trignometric functions as suggested 
by considering the similar integral 



I 



•2/ \ 

sm [nx) 
dx t. — = n 



dy 



sin 2 (y) 



, sin 2 (i/) ^, 
dy £1 = 0(n +1 ) 

y 



(61) 



To find the scaling, we have summed E n over n up to a 
cutoff N using 



N 



sin 2 (no;) sin 2 (ra/) 



M sin(Ma;) sin(My) 
8 8 sin x 8 sin y 

sin(Mx-) sm(Mx+) 
16sinx_ 16sina;_|_ 



(62) 



where M = 2N + 1 and x± = x ± y. Then the remain- 
ing integrations arc done numerically on Mathematica to 
evaluate 



n=l 

In Fig. [3 we plot P{N) versus N for the loop with A = x 
and B = y, clearly showing that the total power radiated 
grows lineary with the upper cutoff N. 

Since P(N) oc N, the total power emitted diverges 
linearly and 



P(N) ~ j^N 



(63) 



4000 



3000- 




1000 



FIG. 7: The plot of 2it{-k 2 L/ Ae) 2 P(N) versus N for the de- 
generate kinky loop shows P(N) oc N growth. 



The kink width provides a cutoff on harmonics and we 
assume this is similar to the width of the string, 1 / ^/TJ. 
Therefore it is reasonable to take N ~ \JJiP and estimate 
the total power emitted, 



E 



2 ( 1 

e n 



(64) 



This energy emission rate is much larger than the naive 
dimensional analysis result e 2 /L 2 , by a factor of ^[JiL ~ 
L/w where w is the width of the string. 



B. Cuspy Loop 

We now consider smooth loops that have cusps i.e. 
points that reach the speed of light . 
A simple loop trajectory is given by 

-(sin s_, — cos s_, 0) 



2tt 



6° 



(T+, b 



L 
27 



(sin s+, 0, — cos s+) 



(65) 



where s± = 2ira±/L. The loop has cusps because there 
are points such that a' = —b'. These occur at s + = s_ = 

0. 7T, at which point the velocity of the string is v = ±x 

1. e. the point on the string reaches the speed of light. 

Inserting the loop trajectory in Eq. (|49p gives 



r(n) 



(-l)V n ^-L 



J n {jlKxy) 



j(n) = e -tn<t, + L 



iJ' n {nK xz )k 



(66) 



where k^ 1 = Lk 1 * / (&irri) , 

K v7 = (k x ,0, Kz) , Kij = (k^jK^jO) 



«4 = ("«*> °! «*) : «x U = -«i> 0) (67) 
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and we have defined the unit vectors k xz , k xy , and 



tan = — , tan0_ = — 



(68) 



The rate of energy emission is obtained from Eq. (pi]) . 
We rescale all momenta and energy by L/ Ann and obtain 



2 4 />1 

en 



2ttL 2 



dq / d 2 p / d 2 p' q(l - q) X 



(OP - (1 - q)P') ■ ( J " «sy - iJ'n( nK xv) k xT, 



(69) 



where k = qp + (1 — q)p' . 

Note that the integrand is non-singular because of the 
properties of the Bessel functions (see Fig. [2]) . The inte- 
gration is over a compact region and can be evaluated 
numerically. For the evaluation it is advantageous to 
rewrite Eq. (joT))) using the Bessel function identities 



2 n J ri{x) _ j n _ 1 ( x ^ _|_ J n+l (x) 



2J' n (x) = J n -l(x) - J n +l(x) 



(70) 



(71) 



To estimate the dependence of E n on n, we compare 
the present integration to that for the kinky loop. The 
main difference is that the integrand for the kinky loop 
contains trignometric functions whose argument contains 
n, whereas here we have Bessel functions. As larger val- 
ues of n are considered, the Bessel functions contribute 
for smaller windows of n xy and k Kj whereas the peri- 
odic trignometric functions contribute for smaller win- 
dows but there are more such windows because of the 
periodic nature of the trignometric functions. Following 
the argument below Eq. (|60[) . we expect E n oc n° and 
hence P(N) oc N for the cuspy loops. 

We have numerically evaluated the integration in 
Eq. dHll) and, in Fig. [SJ we plot 8irL 2 E n /e 2 versus n = 
1, . . . , 10. The plot confirms the E n oc n° estimate for 
the cuspy loops. 

The string thickness, w ~ provides a natural 

cutoff on the wavelength of the emitted radiation. Thus 
the maximum harmonic ~ L/w and we estimate the total 
AB power for massless particles from the cuspy loop 



E 



L 2 



e 2 L 
L 2 w 



(72) 



This is of the same order of magnitude as the radiation 
from kinks. 

The emission of massive particles is conceptually sim- 
ilar but technically more involved. We have left that 
analysis for future work. 
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FIG. 8: (8nL 2 /e 2 )E n versus n for the loop with cusps. 



Identify 



FIG. 9: A conical space is obtained by cutting out a wedge 
from a Euclidean two dimensional space and identifying the 
edges. The wavefunction of a particle that goes around the 
tip of the cone, acquires a phase given by the deficit angle. 



V. GRAVITATIONAL AB RADIATION 

The above analysis for AB radiation suggests that a 
corresponding effect must exist in a gravitational setting. 
Consider the conical metric of a straight cosmic string 
(see Fig. |9]). When any particle goes around the string 
and returns to its original position, the wavefunction ac- 
quires a phase factor given by the conical deficit angle 
5 = 8ttGu, where u is the string tension. Hence any 
particle will interact with a cosmic string by the grav- 
itational AB effect and we expect an oscillating cosmic 
string to radiate gravitational AB radiation that contains 
all species of particles. 

The existence of gravitational AB radiation may also 
be argued on the grounds that an oscillating cosmic 
string provides a time dependent background for quan- 
tum fields. All the vacuum modes of the quantum fields 
have to adjust to the dynamical background. Then it is 
natural to expect some "sloshing" i. e. pair creation. The 
process has been studied in 0] and here we extend that 
analysis to the case of kinks. 

The most intriguing aspect of gravitational AB radi- 
ation is that it implies that cosmic strings will radiate 
photons. Based on our calculations of the previous sec- 
tions, we expect the emission to be largest from kinks on 
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loops. We will now sketch a calculation of the gravita- 
tional AB radiation of massless scalars from a degenerate, 
kinky loop. The relevant interaction between the metric 
and the scalar field can be obtained if we start with the 
action for a massless scalar field 



S,l, 



and expand around the Minkowski metric, 

9^ = V„u + v . 9^ = vr - 

This gives 



Sit, — 5*0.0 



d^x^H^d, 



(73) 



(74) 



(75) 



where S^o is the action in a Minkowski background, and 




1000 



(76) 



The metric perturbation, /i M „, is sourced by the string 
and can be written in Fourier space as 



IQnG 



nh (») 



1 



ji(s)a 



(77) 



where the superscript (s) on T denotes that this is the 
Fourier transform of the string energy-momentum tensor 
which is given by [l2[ 

T$(x) =fij dPa{X^X v - X^XDSHx - X) (78) 

Now the string energy-momentum tensor can be written 
in terms of a(cr_) and 6(<7+). 

rW"" = -~ J d 2 o±{a> l 'b v ' + a v 'b^)5 4 {x - X) (79) 

The Fourier transform can be done and it is clear that it 
will factorize into integrations over left- moving (cr_ ) and 
right- moving (o+) variables. For our present analysis, 
we would like to focus on the degenerate, kinky loops, 
for which a^' and b^ are piecewise continuous. Then 



where I± have the same meaning as in Eq. (|43[) . 
In momentum space we have 



(80) 



WttG - 

U2 r v 



(81) 



and the pair production amplitude from the square loop 
is 

iMW = -^TW(pV'+?y') 



I . p 



7 (p-I+p' -I-+P-I-P' ■/+) (82) 



FIG. 10: Plot of the power radiated, 27r(7r 2 L/85) 2 P(iV) up 
to some harmonic N for gravitational AB radiation from the 
square loop, showing the P(N) oc N growth. 



We use the orthogonality relations in Eqs. f|44[) . together 
with energy conservation, to simplify the amplitude ex- 
pression as in Sec. IIV A| to eventually obtain 



\M 



4^1-1 q\l-qf[q-Aq-B] 



sin 2 ((1 + k ■ AY-f) sin 2 ((!-«• B)f) 



[1 - (k • A) 2 ] 2 [l - (#s • B) 2 ] 2 
where, as in Sec. IIV A\ 

L(p + p') 



K = 



Ann 



and 



Q = 



M 

Ann 



<1 



P ~P 

W - p\ 



(83) 



(84) 



(85) 



The expression in Eq. (|83p for gravitational AB radia- 
tion from square loops is very similar to the expression for 
gauge AB radiation from degenerate, kinky loops, given 
in Eq. (|5T[) . The main difference is in the terms preced- 
ing the trignometric functions. The [q ■ A q ■ B} 2 factor 
implies that gravitational AB radiation is in the plane of 
the loop, in contrast to gauge AB radiation. However, 
the scaling with n follows the discussion below Eq. (|60p , 
and is expected to be independent of the form of the pref- 
actor. As a check on the argument, we have evaluated 
P(N) numerically and the results are shown in Fig. [TOl 
confirming that P(N) oc N i.e. E n oc n°. 

So far we have outlined gravitational AB radiation into 
massless scalars. We can proceed similarly for the radia- 
tion of photons. Now the linear order interaction is 



d 4 x 







(86) 
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where the field strength is defined in terms of the gauge 
potential as = d^A v — d y A^ and indices are raised 
and lowered by the Minkowski metric. The rest of the cal- 
culation is expected to be very similar to that for massless 
scalars, though the final state will also depend on pho- 
ton polarizations. We expect that there will be string-like 
signatures contained in the polarizations, e of the emit- 
ted photons, e.g. linear polarization along the string due 
to factors like e ■ A in the amplitude, and these would be 
worth exploring further. 

As in the gauge AB radiation case, the P(N) cx N (i.e. 
E n (x n°) result suggests that the power radiated in pho- 
tons due to the gravitational AB effect can be estimated 
as 



Eab 7 ~ 5 2 jU 



1 



(87) 



The loop also emits power in gravitational waves, 



E 



Hence 



gw 



Eab 7 ~ E g 



5 fx 



(89) 



where w denotes the width of the string. Even though 
the power emitted in photons is a very small fraction 
of the power emitted in gravitational waves, it can still 
be of interest, especially because it is in photons. For 
G/i = 10~ 8 , we can write /i ~ 10 51 ergs/s, from which 
we get Eg ~ 10 43 ergs/s and E 1 ~ 10 35 (ui/L) ergs/s. 
The lifetime of the loop is determined by its gravitational 
radiation rate and is estimated as ~ L/S. Hence the total 
energy emitted in photons is ~ With 5 ~ 10~ 8 , 

the total energy emitted in photons ~ 10 7 GeV, and is 
independent of the size of the loop. 

The number of photons emitted per unit time in the 
n th harmonic, N n , is proportional to n _1 . Hence the 
same number of photons are emitted in every logarithmic 
frequency interval 



62 

dN u ~ — — 
L oj 



(90) 



In [8[ the accummulated photon background due to 
cusps on a network of cosmic strings was considered. 
Here we estimate photon emission from kinks along a 
single long string in the present universe. This is rele- 
vant for direct detection of a string. The signal may be 
stronger than that from cusps particularly because kinks 
are expected to be much more numerous on strings as 
compared to cusps. A recent estimate of the number 



density of kinks on a string is [14 
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t W* 



(91) 



where t* is the epoch at which friction on the string net- 
work dynamics ceases to be important compared to the 



Hubble damping. Based on the scattering rate of parti- 
cles off the string, a rough estimate is ~ tf/S where t f 
is the time of formation of the strings. 

To obtain numerical estimates, we have the present 
cosmological epoch at to ~ 10 17 s, and for strings with 
5 = 1CT 8 , we find t* - 1CT 27 s. Thus t /U ~ 10 44 . 
With each kink emitting a total energy of over a 
Hubble time, the emission per unit length of string per 
unit time is ~ SyfjKto/t*)/^ ~ 10 4 crgs/cm-s, which is 
~ 10 25 ergs/kpe-s. The energy emitted from strings on 
astrophysical length scales is quite small compared to as- 
trophysical sources (e.g. the solar luminosity is ~ 10 33 
ergs/s), though the unique stringy features, (such as lin- 
ear emission, polarization, flat spectrum) of the emission 
may help observations. Even if the signal is too weak 
to be used to find strings on the sky, it might be a use- 
ful signature to confirm the presence of a string if one is 
suspected in some field of view. 

Before closing this section it is important to point 
out that our estimate treats all kinks equally whereas 
Copcland and Kibble [l4| emphasize that most kinks are 
nearly straight and it is relevant to consider the sharp- 
ness of the kink. Hence a more rigorous estimate should 
fold in the kink sharpness distribution with the photon 
emission rate, and also include effects of kink smoothing. 
(The analysis in 14j accounts for string straightening due 
to Hubble expansion but not kink smoothing due to radi- 
ation backreaction.) We leave more realistic calculations 
of observational signatures for future work. 



VI. CONCLUSIONS 

We have studied the problem of an oscillating solenoid 
in vacuum and determined the rate of scalar particle 
creation due to the Aharonov-Bohm interaction. This 
is novel because the electromagnetic fields of a moving 
solenoid vanish everywhere except at the location of the 
solenoid. In particular, a moving solenoid does not emit 
electromagnetic radiation but does emit charged particle 
pairs as AB radiation. We have also determined the pat- 
tern of AB radiation and find that it is pre-dominantly 
in the direction perpendicular to the plane of oscillation. 

Cosmic string loops that have AB interactions with 
quantum fields will also emit AB radiation. In the thin 
string limit, loops with cusps and kinks emit a linearly di- 
vergent power. Imposing a cutoff due to the thickness of 
the string indicates that the AB radiation is significantly 
stronger than the naive dimensional estimate. 

We point out that all quantum fields interact with cos- 
mic strings via the gravitational AB effect, and thus there 
is corresponding gravitational AB radiation of all parti- 
cles, including photons. Cosmic strings, even if occurring 
in dark matter sectors of particle physics, will emit pho- 
tons. We have calculated the energy emitted in photons 
from loops with kinks and cusps. As in the gauge case, 
the radiated energy is much larger than the naive di- 
mensional estimate. We have made a rough estimate of 
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the energy emitted in photons and it remains to be seen 
if the signal is strong enough to make it an interesting 
observational tool. 
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Appendix A: Moving Frame Perturbation Theory 

Consider a system of coupled oscillators governed by 
the time-dependent Lagrangian 

L = \Y,( d ^) -\Y, K ^ x ^ ( A1 ) 

i ^ ' ij 

The time-dependence enters through the coupling matrix 
K{t). It is helpful to determine the instantaneous normal 
modes q a {t) with instantaneous frequencies uj a (t) that 
satisfy 



(A2) 



These modes are normalized according to 

y^gmgffi = 8 a p- (A3) 

Z 

The oscillator co-ordinates may then be expanded 

Xi = Y,Za{t)q ai {t) (A4) 



where £ a (t) are instantaneous "normal mode" co- 
ordinates. The Lagrangian cq (|A1I) may be written in 
normal mode co-ordinates as 

^E(f (f + -w,)-iv>^. 

(A5) 

where a summation over the repeated indices a and 7 is 
left implicit. Here the "transition element" 



(A6) 



is determined by the evolution of the instantaneous nor- 
mal modes. Eq. (|A5|) generalizes the textbook treatment 
of normal modes for coupled oscillators to the case that 
the Lagrangian is time-dependent. In the limit that the 
coupling matrix K is time-independent, A = 0, and eq 
(|A5[) reduces to the textbook normal mode Lagrangian. 

From the normal mode Lagrangian eq (|A5[) we may 
pass to the Hamiltonian 



H 



2 ^ 



IT 



2 t 21 



a/3 



A a p£,p 



where the canonical momentum is 
Hamilton's equations then read 



dC 



(A7) 



(A8) 



— j = n Q - A a p£p, 



df 

dIL a 
dt 



J2u fj Afs a -u>lt, a . 



(A9) 



In the limit that the coupling matrix K is time- 
independent, A = 0, and the normal mode co-ordinates 
simply undergo the expected harmonic oscillations. In 
the limit that the coupling matrix varies "slowly" , the 
normal modes also vary slowly, allowing us to solve 
Hamilton's equations by treating A as a perturbation. 
To zeroth order the solution is 

£i 0) W = U-T)cos[cj a (t + T)} 
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Ha(-r) 



aa[u a (Jt + T)], 



n« 0) (i) = n Q (-T)cos[^ Q (t + r)] 

- u a ia(-T)mD[u a (t + T)]. (A10) 



The first-order correction is 



n«(t) 



dr ^ ; (t)A^ q (t) cos[w Q (t - t)] + 



/ 1 
dr— IE (r)A /3Q (r) sin[w Q (t - r)] 
-T u o 



dr nf{T)Af, a {T) cosK(i-r)] 



— T 
i 



dr w a £]p (r)Ap a (T) sin[a; Q (t - r)] 



(AH) 



For simplicity in the perturbative solution Eqs. (|A10|) 
and (|A11[) and hereafter it is assumed the eigenfrequen- 
cies u) a (t) are independent of time although the eigen- 
modes q a (t) may be time-dependent. It is not difficult to 
obtain the solution in the case of time-dependent eigen- 
frequencies but the formulae are more cumbersome and 
for the application to solenoids the special formulae given 
here suffice. This is because as the solenoid moves the 
eigenmodes shift but their frequencies do not change. 

Thus far we have focussed on the solution to the clas- 
sical equations of motion. However for a linear system 
such as Eq. (|A1[) the same solution applies in quantum 
mechanics if we interpret n a (t) and £ a (t) as operators in 
the Hciscnbcrg picture that are solutions to the Heisen- 
berg equations of motion. 

To further analyze the quantum mechanics we intro- 
duce the ladder operators 



aUt) 



/ 2uj a 
i 



U a {t) 

n Q (*). (A12) 



Eqs. (|A12I) . (IA10|) and (|A11|) show that the creation and 
annihilation operators at time t are related to those at 
the earlier time —T via the canonical transformation 

a a (t) = ^ up a ap(-T) + vp a aJp(-T), 



(A13) 



where the Bogolyubov coefficients 



up 



x dr Ap a (r)e 

J-T 



x / dr Ap a (r)e- i ^ +u ^ T . 



(A14) 



Eqs. (|A13[) and (|A14[) contain complete information 
about quantum particle production via the dynamical 
Aharonov-Bohm effect. For example the number of par- 
ticles in mode a at time T, (n a (T)) — (a Jf a (T)a a (T)) is 
given by 



-T 



dr Ap a {T)e- %{uJa+u ? )T 



If we take the transition matrix element to oscillate in 
time with frequency f2, 



Ap a {r) = ap a cosflr 
then the rate of particle production in mode a 



(A15) 



Ro 



-y 
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| ap a | 2 S[Q - (u a +wp)]. 

(A16) 

Here we have made use of the usual golden rule prescrip- 
tion T = nS[fl — (u a + cop)]. 

Now let us adapt this general formalism to the case 
of a moving solenoid. In place of Eq. (|A1[) we have the 
Lagrangian 



C 







(V - ieA)0 



(A17) 



that describes a massive scalar field coupled to the vector 
potential of the solenoid, which is given by Eq. (|2"5j) in 
the singular gauge we use here. In place of Eq. (|A2[) the 
instantaneous eigenmodes satisfy 

- [V - ieA{r, t)] 2 ^ a (r, t) + m 2 # Q (r, t) = w 2 * Q (r, t). 

(A18) 

If the solenoid is located along the z-axis the eigenmodes 
are given by Eq. p0[) . Here we consider a solenoid mov- 
ing in the x-z plane while remaining parallel to the z- 
axis; at time t its x-displacement is (vq/Q) s'mQt. The 
eigenmodes at time t, ^ a (r,t) are therefore obtained by 
simply displacing the eigenfunctions Eq. ([30jl along the 
x-axis by the appropriate amount. Thus the modes are 
labelled by a — > k, q, I and have eigenfrequency 



Wfc,,,l(*) = Vk 2 + q 2 +' 



(A19) 



Note that for this problem the instantaneous mode fre- 
quencies do not vary in time although the mode functions 
do vary in time. 

Now the transition elements Eq. (|A6[) are given by 



Ak,q,l;k',q',l' (t) 



f d 

v cosflt / dr V fcj9j i(r)^Vfe',g',i'(r), 
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leading to Eq. (|35|) for the reduced transition element 
a k,q.i-,k' .q' ,i' given in the main body of the paper. The ex- 
plicit form of ak,q,i;k' ,q' ,l' worked out there [Eqs. (|36|) and 
([57)) ] reveals that particles are produced only in modes 
with I = or I = 1. 

Thus far we have focussed on classical analysis. Since 
this is a complex field to analyze the quantum field theory 
we must introduce two sets of creation and annihilation 
operators per mode: c\, t and c k ,q,i that create and an- 
nihilate particles in mode k, q, I and d k l and d kq i that 
create and annihilate anti-particles. 

It follows from Eq. (|A16|) that the rate of particle pro- 
duction in mode (k, q, I) is given by 




X | a kq l. lk 'q'l' | 2 S p - (LJ kq l + UJk'q'l')} (A20) 

Recall that we had discretized our modes for convenience 
by placing the scalar field in a box of dimensions a and 
L. Taking the continuum limit a — > oo and L — > oo in 
Eq. (|A20[) . summing over the azimuthal quantum number 
I yields the expression for I(k,q) in Eq. (|55)1 . the rate 
for production of a single particle species in modes with 
wave- vectors k, q. In taking the continuum limit we make 
use of the usual rules 
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